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1 INTRODUCTION 



Gigahertz Peaked Spectrum (GPS) radio sources are a class 
of powerful compact extra-galactic radio source charac- 
terised by an overall turnover in their radio spectra at about 
1 GHz in frequency (for a review on GPS sources, see O'Dea 
1998). They typically exhibit low flux density variability, im- 
plying that Doppler boosting is probably not important in 
these objects. This view is strengthened by the VLBI mor- 
phologies - of at least those identified with galaxies (not 
GPS quasars) - which tend to show symmetric structures 
on both sides of a central core. These sources are therefore 
also called Compact Symmetric Objects (CSO). The overlap 
between the two classes is large, but may not be complete 
(see Snellen et al. 2000a for a discussion). 

The characteristics of GPS sources are quite different 
from that of other classes of compact radio source, which 
can be highly variable, and are almost always dominated 
by Doppler-boosted, one-sided core-jet emission. This may 
also be the case for GPS sources optically identified with 
quasars which also show core-jet VLBI morphologies. It has 
been argued that the GPS galaxies and quasars may not be 
related at all, but just happen to have similar radio spectra 
(Snellen et al. 1999). In contrast to GPS galaxies, the GPS 
quasars are almost exclusively found at very high redshifts 
(O'Dea et al. 1991; Stanghellini et al. 1998; de Vries et al. 
1997). 

CSO and GPS galaxies receive considerable attention, 
because they are generally thought to be the young progen- 
itors of large-size radio sources. Measurement of hot-spot 
advance speeds (Owsianik & Conway 1998; Owsianik et al. 
1998; Taylor et al. 2000; Tschager et al. 2000) by multi-epoch 
VLBI observations have proven for a handful of archetype 
GPS sources that these objects are indeed young, typically 
10 2-3 yrs. In addition, measurements of the high frequency 
breaks in the slightly larger Compact Steep Spectrum (CSS) 
radio sources, indicate that they have radiative ages in the 
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range of 10 3-5 yrs (Murgia et al. 1999). It is therefore gen- 
erally believed that GPS galaxies evolve into CSS and then 
FRI and/or FRII radio sources, assuming that they are fu- 
elled for long enough. This makes CSO/GPS galaxies the 
key objects to study the early evolution of powerful radio- 
loud AGN (Fanti et al. 1995; Readhead et al. 1996; O'Dea 
& Baum 1997; Snellen et al. 2000b). 



Complete samples of GPS sources are vital to perform 
statistical studies. The complete sample of Stanghellini et al. 
(1998), containing sources with Ssghz > 1 Jy, is the most 
extensively studied. Other samples include the faint GPS 
sources of Snellen et al. (1998), and that with GPS sources 
from the JVAS survey (Marecki et al., 1999) . Since radio 
telescopes are mostly situated in the northern hemisphere, 
radio surveys are also strongly biased towards the north. 
This hampers complete statistical studies of southern GPS 
sources, even at the highest flux density levels. Since the 
space density of GPS sources is so low, full sky coverage 
is important to firmly determine the behaviour of the high 
end of the luminosity function of GPS sources, which is vi- 
tal for constraining radio source evolution models (Snellen et 
al. 2000b). In contrast, the most powerful optical telescopes 
are found in the south (at least those available to European 
scientists) , and therefore it is of great value to have a signifi- 
cant and well defined sample of GPS sources in the southern 
hemisphere. Pioneering work on southern GPS sources has 
been done by de Vries et al. (1995; 2000), especially on ob- 
taining optical identifications and redshifts, and by King et 
al. (1996) involving VLBI observations, but their samples 
are not clearly defined. This paper describes the selection 
of a complete southern/equatorial sample of GPS galaxies 
from the Parkes PKSCAT90 database (Wright & Otrupcek 
1990), and presents the radio spectra and the optical iden- 
tification status in a homogeneous way. Note that the GPS 
sources identified as quasars are not included in the sam- 
ple. Section 2 describes the selection of the sample and its 
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radio properties. Section 3 describes optical CCD imaging 
and spectroscopy performed with the ESO 3.6m telescope. 
In section 4, the results are presented and discussed, in- 
cluding an update on the R-band Hubble diagram for GPS 
galaxies. 



2 A SOUTHERN/EQUATORIAL SAMPLE OF 
GPS GALAXIES 

2.1 Selection of the sample 

The basis for the selection of sources in the sample is 
the Parkes multi-frequency survey (PKSCAT90, Wright & 
Otrupcek 1990), which is a compilation of flux density mea- 
surements at frequencies of 178 MHz (4C), 408 MHz (Mo- 
longlo catalogue), and at 625, 1410, 2700, 5000, and 8400 
MHz taken from the Parkes surveys. Firstly, all sources 
with 52700MHz > 0.5 Jy were selected in the regions with 
-40° < Dec < +15° and galactic latitude, \b\ > 20° (1129 
sources) . This flux density limit and galactic latitude cut-off 
were chosen to be the same as those of the Parkes half- 
Jansky sample of flat spectrum radio sources (Drinkwater 
1997), allowing comparison of source statistics between the 
two samples. Note that most GPS galaxies are not found in 
the flat spectrum sample, since their high frequency spec- 
trum is too steep. Since 1410 MHz Parkes data is only 
sparsely available, the positions of these sources were cross- 
correlated with the NVSS survey (Condon 1998), of which 
all flux density was included within a radius of 6' from 
the Parkes position. In addition, the sources were cross- 
correlated with the Texas 365 MHz survey (Douglas et al 
1996), which covers most of the sky area of the sample. 
Each spectrum was then fitted in an automated way with a 
2nd order polynomial function to find objects with a possi- 
ble turnover. A total of 323 objects were pre-selected all of 
which showed a maximum above 100 MHz in the resulting 
fit. The spectra of these objects were then inspected visually 
to judge the validity of the fit. The resulting 193 objects 
were cross-correlated with all Parkes radio sources in the 
NASA/IPAC Extra-galactic Database (NED) with available 
redshifts, this to filter out any quasars. The NED database 
and the CATS database of radio sources (Verkhodanov et 
al. 1997) were then searched extensively for additional flux 
density points and other information on the remaining 125 
objects. In addition, the NVSS maps of each of the candi- 
dates were checked to see whether or not any extended emis- 
sion was present. The low frequency survey measurements 
were taken with interferometers, and therefore possible ex- 
tended emission may have been missed. This indeed turned 
out to be a major source of contamination for many (47) 
mostly nearby radio galaxies. In other cases, the extra data 
points found in the NED and CATS databases showed that 
the spectra did not turn over at frequencies below 300-400 
MHz, or were subject to strong variability (19). A few can- 
didates turned out to be planetary nebulae (3), and some 
were BL Lacs or quasars without a measured redshift (7). 
The spectra of all sources which were not pre-selected (806 
objects) were also visually inspected to see whether any can- 
didates had been missed. The remaining sample contains 49 
radio sources with peaked spectra, which are optically iden- 
tified with a galaxy or still without an optical identification. 



2.2 The radio properties of the sample 

From the 49 sources in the resulting sample, 12 are also 
present in the brighter samples of GPS sources from 
Stanghellini et al (1998) and de Vries et al. (2000). The prop- 
erties of the objects are given in table [l], with in column 1, 
the J2000 IAU name, in column 2 an alternative name, in 
column 3 the most accurate radio position found in the NED 
database or that from the NVSS survey (on which the opti- 
cal identification process is based), in column 4 the R-band 
magnitude (corrected to Cousins R and galactic extinction 
taken into account), in column 5 the redshift, in column 
6 the flux density at 2700 MHz, in column 7 the turnover 
frequency, in column 8 the flux density at the turnover, in 
column 9 the reference for the magnitude, in column 10 the 
reference for the redshift, and in column 11a possible com- 
ment. The radio spectra are shown in figure |l|. The solid 
lines indicate a 2nd order polynomial fit to the data-points 
to guide the eye, and do not represent physically meaning- 
ful model-fits. The turnover frequencies given in column 6 of 
table ^ are often just estimates, since it was not possible to 
determine them accurately due to the sparse sampling near 
their spectral peak. 



3 OPTICAL ESO 3.6M OBSERVATIONS 

Optical CCD imaging and spectroscopy was performed on 
objects in the sample using the ESO 3.6m telescope at La 
Silla, in Chile, on 27-28 July 2000, using the ESO Faint 
Object Spectrograph and Camera (EFOSC2). The weather 
was photometric throughout both nights. EFOSC2 uses a 
2048x2048 thinned Loral CCD chip with a 0.157 pixel size, 
which was binned by 2x2 for both the imaging and spec- 
troscopy. 

3.1 CCD imaging 

Fourteen objects in the sample were observed in the g, r, 
and/or i filters, which have effective wavelengths of 5169A, 
6814A, and 793lA respectively. For flux density calibra- 
tion standard star fields from Landolt (1992) were observed 
throughout the nights, and fluctuations in the magnitude 
scale zero-points were found to be smaller than 0.03 mag. 
The images were reduced in a standard way in the software 
package IRAF from the NOAO, but were positionally cali- 
brated within the AIPS software package (from the NRAO) 
using the task XTRAN. Positions of bright stars in the CCD 
fields were taken from the APM catalogue (Irwin et al. 1994) 
or the COSMOS UKST catalogue (Hambly et al. 2001), and 
in some cases directly measured on the 2nd Digitized Sky 
Survey. The flux calibrator measurements in Landolt (1992) 
were performed in the Johnson-Kron-Cousins UBVRI pho- 
tometric system, which is a somewhat different filter system 
from that used in this paper. The differences for the R and I 
band observations were found to be within the observational 
uncertainties, however for the g band observations, a colour 
correction to V-band had to be applied. From the observa- 
tions of the standard star fields, this colour correction was 
found to be V — g— 0.4A(V— R). The magnitudes of the opti- 
cal identifications were determined by adding the flux within 
a box including all the emission from the object nearest the 
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Table 1. The radio and optical properties of objects in the southern/equatorial sample of GPS galaxies. 
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Comments: 

A) Also in the sample of >Uy GPS sources from Stanghellini et al. (1998); B) Also in the sample of de Vries et al. (1997) 
C) Also in the sample of de Vries et al. (2000); D) Possibly a quasar; E) Large offset between the radio and optical position 
F) Near a bright star, no magnitude 
References: 

1) This paper; 2) Snellen et al. (1996); 3) O'Dea et al. (1990); 4) de Vries et al. (1995); 5) de Vries et al. (2000) 
6) Biretta et al. (1985); 7) Stern et al. (1999); 8) Fugmann et al. (1988); 9) Stanghellini et al. (1993) 
10) Stanghellini ct al. (1998); 11) Wright & Otrupcck (1990); 12) Drinkwater et al. (1995) 

13) Digitized Sky Survey II; APM catalogue (Irwin et al. 1994); SuperCosmos Sky Surveys (Hambly ct al. 2001) 
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J0022+0014 J0108-1201 J0206-3024 J0210+0419 J0210-2213 J0242-2132 J0323+0534 




J0401-2921 J0407-3924 J0407-2757 J0433-0229 J0441-3340 J0457-0848 J0503+0203 




J1 122-2742 J1 135-0021 J1203 + 0414 J1345-3015 J 1 347+1 21 7 J1350-2204 J1352 + 0232 




J 1 352+1 1 07 J1447-3409 J 1 506-09 1 9 J1521+0430 J1543-0757 J1546 + 0026 J 1 548-1 2 1 3 




J1556-0622 J1600-0037 J1604-2223 J1640+1220 J1648 + 0242 J1734 + 0926 J2011-0644 




J2058+0542 J2 1 23-01 1 2 J213Q + 0502 J215 | 1 +0552 J22|2 + 0152 J2325-0344 J2339-0604 



Figure 1. The radio spectra of the GPS galaxies in our sample. The tick-marks indicate 0.5, 1.0 and 5.0 GHz in frequency, and 0.5,1.0, 
and 5.0 Jy in flux density. 
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Table 2. Details of the CCD imaging observations. 



Name 




V band 


R band 




I band 




t 


my 


t 




t 


771/ 




sec 


mag 


sec 


mag 


sec 


mag 


J0108-1201 


60 


23.50±0.23 


120 


22.45±0.06 


300 


22.15±0.15 


J0210-2213 


60 


>24.40 


1500 


23.55±0.13 


300 


23.68±0.32 


J0210+0419 






1500 


>24.20 






J0407-2757 


60 


21.94±0.19 


60 


21.24±0.04 






GPS- west 




22.61±0.16 




21.71±0.11 






east 




22.87±0.20 




22.12dz0.13 






J1350-2204 


60 


21.99±0.07 


300 


21.13±0.05 


300 


20.30±0.12 


GPS-north 




23.33±0.29 




22.04±0.08 




20.89±0.08 


south 




22.59±0.12 




21.93±0.10 




21.50+0.13 


J1548-1213 


60 


22.93±0.17 


300 


22.52±0.13 


300 


21.55+0.11 


J1556-0622 


60 


>23.40 


1500 


22.82±0.13 


1500 


22.33±0.06 


OclSt 








23.40±0.15 




22.93+0.13 


GPS- west 








23.80±0.17 




23.26+0.17 


J1604-2223 


60 


21.24±0.05 






300 


18.53±0.03 


J1640+1220 






60 


21.48±0.20 






J2011-0644 


60 


22.81±0.08 


1500 


21.52±0.04 


300 


20.48±0.06 


north 








22.06±0.04 




21.28±0.08 


GPS-south 








22.48±0.07 




20.99±0.07 


J2130+0502 






1500 


22.38±0.07 


750 


21.05±0.03 


J2339-0604 






210 


23.00±0.20 


600 


22.12±0.09 



radio position. In the case of multiple objects near the radio 
positions (possible merging or interacting galaxy systems), 
the magnitudes were calculated for the total system and 
each component individually. 

3.2 CCD spectroscopy 

We used the EFOSC2 grism#06, which has 600 gratings per 
mm, resulting in a dispersion of 2.06 A/pixel and a wave- 
lengths coverage of 3860-8070 A. A slit-width of 1" was used 
resulting in a resolution of 13 A. Usually the slit was ori- 
ented near the paralactic angle. The reduction of the spec- 
tra was carried out in a standard way using NO AO's IRAF 
reduction software. Wavelength calibration was carried out 
using a Helium- Argon lamp. One dimensional spectra were 
extracted by summing in the spatial direction over an aper- 
ture as large as the spatial extent of the continuum or the 
brightest emission line. 

3.3 Digitized Sky Surveys 

If the object was not observed by us, and no optical observa- 
tion was found in the literature, an image was extracted from 
the red 2nd Digitized Sky Survey. In addition, the APM cat- 
alogue (Irwin et al. 1992) and the SuperCosmos Sky Survey 
catalogues (Hambly et al. 2001) were searched for a possible 
identification and subsequently magnitude estimate. Unfor- 
tunately, the galaxy magnitudes in these catalogues are only 
accurate to about 0.5 mag. DSS2 data were retrieved for 22 
sources (see figra), with 9 objects remaining unidentified. 



4 RESULTS AND DISCUSSION 

The photometric results of the imaging part of the EFOSC2 
observations are given in table 0, with in column 1 the source 
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Table 3. Details of the spectroscopic observations. 



Name 


t 


Spectral feature 


z 




sec 


name 






J1548-1213 


1800 






0.883zb0.001 






[OII]3727A 


7019 


0.8833 






[NeIII]3869A 


7285 


0.8829 


J1604-2223 


2700 






0.141zb0.001 






G4300A 


4909 


0.1416 






mgb5169A 


5898 


0.1410 






Na5899A 


6727 


0.1404 






H/34861A 


5549 


0.1415 


J 1640+ 1220 


4800 






1.150zb0.003 






Mgll 2798A 


6015 


1.1497 






CII]2326A 


5010 


1.1539 






CIII]410lA 


1909 


1.1482 


J2011-0644 


4200 






0.547+0.001 






[OII]3727A 


5766 


0.5471 






break4000A 


6150 


0.54 



name, in column 2 and 3 the exposure time and magnitude 
of the V band observations, in columns 4 and 5 those of the 
observations in R band, and in column 6 and 7 those of the 
I band observations. The R band images are shown in fig- 
ure 0, except for J1604-2223, for which the I band image is 
shown. Four objects (J0407-2757, J1350-2204, J1556-0622, 
J201 1-0644) seem to have a close companion object. For 
these objects, the total magnitude of the system is given, 
and those of the individual components, with that closest to 
the radio position indicated with 'GPS' in column 1 of table 
|j| Note that the magnitudes in table ^ have not yet been 
corrected for galactic extinction. 

The results of the spectroscopic part of the EFOSC2 ob- 
servations are given in table ^, with in column 1 the source 
name, in column 2 the exposure time, and in column 3, 4, 
and 5 the name, observed wavelength and derived redshift of 
the identified spectral features. Only the galaxies for which 
spectral features could be identified and subsequently their 
redshift determined, are listed. Their spectra are shown in 
figure ^. Spectra of several other galaxies were taken, but 
did not show any identifiable features, and therefore no red- 
shift could be determined. Clearly these objects do not have 
emission lines of high equivalent width. 

The ESO results and those of the DSS2 are incorpo- 
rated in table [j], with in column 4 the R band magnitude, 
corrected for galactic extinction using the values of Schlegel 
et al. (1998). The references to the optical magnitude and 
the redshift are given in column 9 and 10. The last column of 
table [l] gives possible comments. One source, J 1109+ 1043, 
has a large offset between the radio and optical position, 
casting doubt on its possible optical identification. The radio 
source J1600-0037 is located too close (17") to a 12th magni- 
tude star to allow an identification. This object is therefore 
excluded from statistical analysis. The optical identification 
of J1135-0021 could be point-like and therefore may be a 
quasar. Strictly speaking it is not clear whether the empty 
fields are actually galaxies and not quasars. However, from 
the magnitude distribution of other samples of GPS sources, 
it is rather unlikely that they are quasars (most of the GPS 
quasars have ibr < 21). In general, only when all the sources 
are identified and spectra taken, can we be confident that 
all the objects are galaxies. The total identification fraction 
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J1548-1213 z = 0.883 J1604-2223 z = 0.141 
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J1640+1220 z=1.150 J2011-0644 z=0.547 
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Figure 2. Optical spectra taken with EFOSC2 on the ES03.6m telescope. Identified spectral features are indicated 



of the sample is ~80%. About half of the identified objects 
have their redshifts determined. 



4.1 The magnitude distribution 

The apparent magnitude of a (GPS) galaxy is closely related 
to its redshift, as will be discussed later. Therefore, the mag- 
nitude distribution of a GPS galaxy sample is a reflection of 
its redshift distribution. This distribution is very interest- 
ing, since it is closely related to the birth-function of radio 
galaxies and the cosmological evolution of the GPS lumi- 
nosity function (Snellen et al 2000). The cumulative mag- 
nitude distribution is shown in figure |^, with the shaded 
area indicating the distribution of the half-Jansky sam- 
ple (taking into account the uncertainty of the magnitudes 
of the objects without optical identification). GPS sources 
with S2.7GHZ > 1 Jy were excluded to avoid overlap with 
the bright sample. The dashed line shows the magnitude 
distribution of the GPS galaxies of the >Uy (at 5 GHz) 
Stanghellini et al. (1998) sample, which is completely iden- 
tified, and the dotted lines indicate the limits of the mag- 
nitude distribution of the galaxies in the faint GPS sample 
of Snellen et al. (1998). The mean flux density of the GPS 
sample of Stanghellini is about a factor of 4 brighter, and 
that of the faint GPS sample a factor of 6 fainter than the 
half-Jansky sample (with S2.7GHZ < 1 Jy). Figure ^ clearly 
shows that radio samples containing gernerally fainter ra- 



dio sources are associated with galaxy samples with fainter 
magnitudes, with the median magnitudes of the bright, half- 
Jansky, and faint samples being about 20, 21, and 22 in R 
respectively. This seems to indicate that the radio-fainter 
GPS galaxies are primarily at higher redshift, and not sim- 
ply lower in their radio luminosity. 

4.2 The Hubble diagram, optical colours and 
absolute magnitudes 

The Hubble diagram of GPS galaxies has been previously 
discussed by Snellen et al. (1996a, 1996b) and O'Dea et al. 
(1996). Our group has claimed that the apparent magnitude 
- redshift relation of GPS galaxies show a low dispersion, 
and that they are apparently fainter at high redshift than 
3C objects. This has been interpreted in terms of the GPS 
galaxies being intrinsically redder than 3C galaxies, possibly 
due to the lack of the extra, blue light associated with the 
radio-optical alignment effect (McCarthy et al. 1987; Cham- 
bers et al. 1987; Best et al. 1997) . This view is strengthened 
by observations that the optical to near-infrared colours of 
GPS galaxies are consistent with the light of these hosts 
being dominated by old stellar populations (Snellen et al. 
1996b; O'Dea et al. 1996). In addition, Best et al. (1998) 
show that 3C galaxies are typically up to a couple of magni- 
tudes bluer in R-K magnitudes than expected for passively 
evolving elliptical models. The lack of luminosity evolution 
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Cumulative magnitude distribution 




16 18 20 22 24 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

R band magnitude Redshift 



Figure 3. The cumulative R band magnitude distribution of the 
half-Jansky sample of GPS galaxies, with the grey area indicat- 
ing the uncertainty caused by the empty fields. The dashed line 
indicates the same distribution for the bright Uy Stanghellini et 
al. (1998) sample, and the dotted lines the limits to that of the 
faint GPS sample of Snellen et al. (1998). 
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Figure 4. The Hubble diagram of GPS galaxies (solid squares) 
and 3C galaxies (open circles). The dotted line indicates the Hub- 
ble relation for GPS galaxies as found by Snellen et al. (1996). 
The crossed GPS galaxies are those new to the sample. 



Figure 5. The R-I colours as function of redshift for GPS galaxies 
taken from Stanghellini et al. (1993), Snellen et al. (1998) and this 
paper. 



found out to z=l, which would be expected for passively 
evolving ellipticals, was interpreted as a possible dynami- 
cal evolution of the systems through accretion, as found for 
brightest cluster galaxies (Aragon-Salamanca et al. 1998). 
Note that most faint (R~24) GPS galaxies are found to have 
R-K~5-6 (Snellen et al. 1996b, 1998), and can therefore be 
classified as extremely red objects (ERO). Given their rela- 
tive ease of selection, GPS galaxies promise to be excellent 
objects for studying the evolution of massive ellipticals out 
to z > 1. 

The combined samples of Stanghellini et al (1998), 
Snellen et al. (1998), and that presented in this paper pro- 
vide an additional 15 GPS galaxies with available redshifts, 
compared to the original analysis of the R band Hubble dia- 
gram of Snellen et al. (1996). A few galaxies have only rough 
estimates of their R magnitude from the DSS surveys and 
are excluded from the analysis below, leaving 34 galaxies in 
total (the galaxy J2151+0552 was excluded as well, since 
its optical light is strongly influenced by non-thermal emis- 
sion). The updated Hubble diagram of GPS galaxies (solid 
squares) is shown in figure ^, in which all magnitudes are 
converted to Cousins R. The dotted line indicates the Hub- 
ble relation as fitted to the original dataset in Snellen et al. 
(1996). The open circles indicate 3C galaxies from a com- 
pilation of samples from Best et al. (1997), Eales (1985), 
and de Koff et al. (1996) with, if necessary, the magnitudes 
converted to 'total magnitudes' in R Cousins. As found be- 
fore, the distant GPS galaxies are apparently fainter than 
3C galaxies at similar redshift, which can be accounted for 
by the extra light of the alignment effect in 3C galaxies (eg. 
Best et al. 1997). However, some differences from the original 
picture now emerge. On average, the nearby objects seem to 
be fainter than in the original sample, and the galaxies at 
z > 1 seem to be brighter than the original Hubble diagram 
indicates. 
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Figure ^ shows the R-I colours as function of redshift, 
with sources taken from Stanghellini et al (1993), Snellen 
et al. (1999) and this paper (all corrected to Cousins R). 
The shaded area indicates the expected R-I colours for a 
passively evolving elliptical with a formation redshift of z=5, 
and a range from 1 to 2 times solar metallicity (Worthey et 
al. 1994). The R-I colours are more or less as expected for 
passively evolving ellipticals, as shown before by O'Dea et 
al. (1996). 

Figure ^| presents the absolute magnitudes for GPS 
galaxies (solid squares) and 3C galaxies (open circles). A 
cosmology with Ho— 70 km sec _1 Mpc _1 , SI a — 0.3 and 
fl m = 0.7 is assumed, and the magnitudes are K-corrected 
using a synthetic galaxy spectrum by Worthey (1994) for 
a 12 Gyr old elliptical with solar metallicity. The models 
of Worthey (1994) are also used to calculate the grey area, 
which represents the expected luminosity evolution as func- 
tion of redshift for a 2L„ elliptical with a formation redshift 
of z=5, and a range of metallicities from 1 to 2 times solar. 

As for the Hubble diagram, this picture seems to be 
different to that found earlier for the absolute K-band mag- 
nitudes of GPS galaxies as function of redshift (Snellen et al. 
1996b), which show a more or less constant absolute mag- 
nitude as function of redshift. We believe this has several 
reasons. Firstly, the few low redshift GPS galaxies avail- 
able in 1996 seem to have been, on average, over-luminous 
with respect to the general low redshift population. This 
was already a concern at that time, since two of them 
have a Seyfert nucleus (4C+12.45 and B1404+286). This 
has moved the z—0 baseline to a lower luminosity, making 
the higher redshift objects relatively brighter. Secondly, the 
now popular A-cosmology makes the objects at redshift of 
unity intrinsically about half a magnitude brighter as well. 
Thirdly, the newly measured redshifts at z > 1 are strongly 
biased towards the objects with the strongest emission lines, 
making contamination from an active nucleus more likely. It 
is therefore expected that the majority of the objects at 
z > 1, for which we were not able to determine a redshift, 
are intrinsically less luminous than the few objects shown 
in figures ^ and []. Indeed, from the four spectra shown in 
figure ^[ the emission lines of the two galaxies at low redshift 
have very low equivalent widths, while those of the two at 
high redshift have very high equivalent width. For that rea- 
son only the redshifts of the latter two could be determined. 
Several deep spectra of other faint galaxies were taken with 
no lines detected. Eight-metre-class telescope time will be 
needed to resolve this issue. 



-22 



-24 



-26 




Redshift 

Figure 6. The absolute magnitudes of the host galaxies of GPS 
sources (solid squares) and 3C galaxies (open circles) as func- 
tion of redshift. The magnitudes are k-corrected using a synthetic 
spectrum of a 12 Gyr old elliptical with solar metallicity from 
Worthey (1994). The crossed GPS galaxies are those new to the 
sample. 



of interaction and/or mergers are larger than that of the 
lifetime of a radio source. 

The statistics in the sample presented here are not yet 
good enough to perform a similar statistical analysis. How- 
ever, it is interesting to note that from the 12 galaxies we 
imaged, 4 are clearly part of a multiple system (J0407- 
3924, J1350-2204, J1556+0622, and J201 1-0644). In addi- 
tion, J1548-1213 seems to have a 2nd component, but the 
picture is less clear due to 2 foreground stars nearby. At first 
glance, this gives a similar fraction of multiple systems to 
O'Dea et al. (1996) and Heckman et al. (1986). It is also 
interesting to note that in all 4 of the multiple systems, the 
GPS radio sources are identified with the redder component 
of the two. A more extensive observing programme for op- 
tical imaging and spectroscopy is planned. 



4.3 Mergers and interactions 

Previous work has indicated that disturbed optical mor- 
phologies and the presence of close companions are com- 
mon characteristics of GPS galaxies (Stanghellini et al. 1993, 
O'Dea et al. 1996). Using a sample of radio- bright GPS 
galaxies, O'Dea et al. find that about half of the galaxies 
have disturbed morphologies and about one-fifth show a 2nd 
component within a few arcseconds. Very similar percent- 
ages have been found by Heckman et al. (1986) for a general 
sample of powerful radio-loud AGN. If GPS sources evolve 
into large scale radio sources, then it would be expected that 
their galactic environments are similar, since the timescales 



5 SUMMARY AND CONCLUSIONS 

In this paper we have described the selection of a new south- 
ern/equatorial sample of Gigahertz Peaked Spectrum (GPS) 
radio galaxies, and presented subsequent optical CCD imag- 
ing and spectroscopic observations using the ES03.6m tele- 
scope. The sample consists of 49 sources with —40° < S < 
+ 15°, \b\ > 20°, and S2.7GHZ > 0.5 Jy, selected from the 
Parkes PKSCAT90 survey. About 80% of the sources are 
now optically identified, and about half of the identifications 
have their redshifts determined. The R-band Hubble dia- 
gram and evolution of the host galaxies of GPS sources are 
reviewed. The view on GPS galaxies seems to be slightly al- 
tered, with the host galaxies increasing in optical luminosity 
with redshift, as expected for passive evolution. This change 
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appears to be caused by the low-z galaxies being less lumi- 
nous than previously thought, the new A-cosmology used, 
and a strong observational bias towards high redshift GPS 
galaxies being those with strong emission lines. 

The sample presented in this paper has the potential 
to greatly improve the statistics on luminous young radio- 
loud AGN, and the use of hosts of GPS sources for galaxy 
evolution studies, especially with the VLT and Gemini-south 
telescopes. 
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A Parkes half-Jansky sample of GPS galaxies 




Figure 8. Images retrieved from the red 2nd Digitized Sky Survey 
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